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1.  Flutter  simulation  is  a  closed  loop  ground  vibration  to  find  flutter  speeds 
on  actual  missile  wings.  The  elements  in  a  closed  loop  convert  the  wing 
linear  and  angular  motions  into  electro-mechanical  forces  which  simulate  the 
aerodynamic  forces  generated  by  the  motions  in  supersonic  flight  at  some 
chosen  Mach  number  and  altitude.  If  the  closed  loop  becomes  unstable,  it 
is  inferred  that  the  missile  surface  would  become  unstable  in  flight. 

RESULTS 


The  results  are  developed  in  terms  of  what  can  be  accomplished  using  one 
shaker  delivering  a  concentrated  force  in  response  to  one  particular  pitching 
angle  on  a  wing.  This  is  a  necessary  first  step  toward  a  system  of  shakers  to 
represent  the  distributed  airforces  which  are  functions  of  the  linear  and  angular 
displacements  and  velocities. 

Studies  of  a  missile  wing  flutter  problem  by  Dr.  E.  Shotland  and  Dr.  A. 
Mitchell  were  made  in  1952  with  the  drastic  simplification  of  the  airforces 
which  is  embodied  in  the  one-channel  simulator  system.  A  check  on  some  two- 
dimensional  translation-pitch  flutter  cases  for  which  classical  NACA  results 
were  available  led  to  the  tentative  conclusion  that  the  lift  force  per  unit 
angle  was  the  most  significant  coefficient  in  the  flutter  problem. 

To  get  a  one-channel  system,  four  major  components  were  needed:  (1)  An 
angular  displacement  transducer,  (2)  a  preamplifier  to  raise  the  signal  level; 

(3)  a  power  amplifier  and  (A)  an  electro-mechanical  shaker.  Various  designs  for 
these  components  were  used  in  the  flutter  simulation  of  two  aluminum  cantilever 
plates,  a  30"  x  12"  x  .125"  plate  A,  and  a  15.5"  x  10"  x  ,25"  plate  B. 

It  should  be  noted  that  component  (1) ,  the  angular  transducer,  posed  a 
durable  problem.  It  was  difficult  on  plate  A  to  match  the  gains  of  two  pick¬ 
ups,  the  differential  output  being  the  pitch  angle  signal,  A  number  of  angular 
transducers  were  tested  t  obtain  quantitative  simulation  for  the  second  plate. 
The  best  device  was  an  electro-optical  angular  pickup  and  an  associated  pre¬ 
amplifier.  For  the  one-channel  system  as  a  whole  a  shaker  force  per  unit  pitch 
angle  up  to  17000  lb/radian  over  a  frequency  band  from  0  to  500  cps  is  now  avail¬ 
able.  It  has  been  used  to  produce  flutter  simulation  of  the  second  plate,  and 
predicts  a  flutter  Mach  No.  of  2.9  at  sea  level.  A  number  of  qualifications 
are  needed  on  this  prediction,  none  of  which  bear  on  the  validity  of  the  one- 
chat  nel  feedback  loop  as  such.  These  qualifications  arise  from  such  issues  as 
the  :hcice  of  suitable  airforce  coefficients  and  the  number  of  channels  which 
are  needed  to  simulate  the  coefficients. 

The  work  serves  to  illustrate  the  performance  of  the  one-channel  system. 

It  is  recommended  that  evaluation  of  the  system  be  continued.  If  it  continues 
to  behave  well,  it  is  recommended  that  another  channel  be  added  for  further 
tests;  a  number  of  channels  should  be  able  to  simulate  both  flutter  and  diver¬ 
gence.  Finally,  it  is  noted  that  the  electro-optical  angular  pickup  should  find 
useful  applications  in  the  measurement  of  static  and  dynamic  angular  motions  of 
structures  in  general. 
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The  topics  which  follow  cover  the  choice  of  airforce  coefficients  for 
simulation,  and  the  details  for  the  two-dimensional  flutter  analysis  using 
the  lift  due  to  angle.  The  modification  due  to  aerodynamic  damping  is  shown. 
Vibration  testing  for  special  angle/force  data  is  described  for  use  in  a 
direct  calculation  of  flutter  speed. 

The  theory  is  extended  to  cover  the  three-dimensional  case  of  a  canti¬ 
lever  plate.  The  concept  for  the  finding  direct  flutter  data  from  open-loop 
vibration  testing  is  extended  to  the  three-dimensional  case.  Closing  this 
loop  then  simulates  the  flutter.  The  development  problems  for  an  early  study 
of  a  Plate  A  and  a  recent  study  of  a  Plate  B  conclude  the  report. 

1.  Choice  of  Airforce  Coefficients  for  Simulation 

The  objective  is  to  make  a  limited  number  of  ground  test  forces 
simulate  the  most  important  aerodynamic  forces.  The  force  which  appears 
to  be  the  most  important  is  the  lift  due  to  angle  of  attack.  This  opinion 
is  reached  by  doing  simplified  studies,  of  the  sort  reported  in  Section  2. 
The  results  are  encouraging.  The  studies  were  originally  made  to  check 
the  assumption  made  by  Dr.  E.  Shotland  and  Dr.  A,  K.  Mitchell  of  the 
Applied  Physics  Laboratory  for  wing  flutter  analysis  in  1951.  Such 
assumptions  allowed  other  complex  factors  to  be  introduced  into  the 
analyses,  and  exposed  the  effects  of  frequency  ratio  and  mass  balance. 
Piston  theory  was  a  later  notable  simplification  by  Ashby  and  Zartarian, 
made  for  the  same  reasons.  In  recent  years,  other  flutter  studies  using 
simplified  airforces  have  been  reported,  and  much  research  is  in  progress 
on  the  subject. 

Tests  at  APL  using  the  one  force  due  to  angle  of  attack  have  been 
made  on  two  cantilever  plates,  A  and  B.  These  tests  prompted  a  review 
of  the  important  airforces,  and  it  is  now  believed  that  damping  forces 
in  the  bending  and  twisting  degrees  of  freedom  are  important.  The  ratio 
of  the  two  generalized  damping  forces  influences  the  flutter  frequency, 
and  hence  the  flutter  speed.  On  plate  B,  (Figures  3  and  4)  the  damping 
force  and  moment  are  provided  by  two  concentrated  forces,  one  3  inches 
ahead  and  the  other  3  inches  behind  the  elastic  axis  or  mid-chord,  1/2 
inch  in  from  the  tip.  The  dampers  are  not  powerful  enough  to  simulate 
the  calculated  generalized  damping  factors  called  for  by  piston  theory, 
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but  they  do  provide  a  step  in  the  right  direction.  The  lift  coefficient 

4  4 

C.  from  piston  theory  ia  —  ,  while  the  classical  value  is  — - 1. 

L  K  V  M2-l 

It  will  be  seen  from  the  discussion,  and  particularly  equation  (32),  that 

either  value  of  C^^may  be  used  to  evaluate  its  effect  on  the  altitude- 

Mach  number  flutter  boundary. 

2.  Flutter  Analysis  -  Two-Dimensional  System 

(a)  One  No-Lag  Lift  Force,  No  Damping 

It  is  well  be  begin  with  a  simple  two-dimensional  surface  having 
bending  and  rotational  degrees  of  freedom,  reacting  only  to  an  aero¬ 
dynamic  force  N a  at  the  elastic  axis,  Figure  1. 

The  equilibrium  equations  for  sinusoidal  oscillation  are: 
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A  condition  of  flutter  exists  when  the  determinant  is  equal  to 
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The  real  and  imaginary  parts  are  separately  set  equal  to  zero: 
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This  value  of  the  frequency  is  then  used  in  the  real  part  of  the 
equation: 
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The  equation  ultimately  simplifies  to: 
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The  equation  can  be  put  into  a  form  suitable  for  correlating  results 
with  NACA  calculations  in  Reference  (1)  by  using  the  definitions  on  page  2. 
and  the  ultimate  result  is: 
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Table  I,  page  27  shows  reasonable  agreement  of  che  results  from 
this  simplified  theory  with  the  "exact”  SACA  calculations. 


(b)  One  lift  force  plus  damping  in  both  degrees  of  freedom 

The  comparison  sho~"a  is  for  a  theory  in  which  the  damping  is 
entirely  structural  of  the  form  i  k  g.  Aerodynamic  damping  is 
viscous  in  nature;  and  an  examination  of  che  flutter  determinant 
shows : 
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The  flutter  frequency  is  determined  from  the  condition  that  the 
imaginary  part  of  the  expanded  determinant  is  equal  to  zero. 


(ID 


(12) 


C 


hT 


h 


+ 


C 


h 


co. 


G  +  C  oo, 
a  a  f 


imm  nmaa  utiunn 


10. 


t 


u 


rJ 

i*} 

if 


The  real  part  of  the  determinant  should  also  vanish: 
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Thus  it  is  possible  to  find  the  critical  value  of  the  aerodynamic 
coefficient  needed  to  satisfy  (13). 
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The  use  of  aerodynamic  viscous  damping  results  in  the  aero-  i 

dynamic  force  coefficient  N  required  for  flutter,  to  be  compared  to 
Equation  (8) ,  where  only  structural  damping  of  the  form  i  k  g  was 
present.  In  any  given  problem,  two  values  of  N  may  be  computed, 
and  the  effect  of  the  aerodynamic  damping  terms  evaluated.  The  major 

f 

effect  arises  from  the  ratio  of  the  two  damping  factors,  to  judge  from  • 

the  example  of  Plate  B.  ! 

The  aerodynamic  coefficients  used  in  the  analysis,  namely  N,  ; 

C.  oo  and  C  Bare  sufficient  to  describe  the  forces  given  by  second  ' 

n  a  os  I 

order  piston  theory  -  Reference  (2).  Then  the  values  are:  > 
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for  forces  at  the  mid-chord  of  a  thin  plate. 


3.  Ground  Vibration  Testing  -  Two  Dimensional  System 

f' 

In  special  cases  it  has  been  shown  that  the  simplified  airforce  theory  y 

with  one  lift  force  yields  reasonable  answers.  Adding  damping  then  in  effect 
makes  it  a  piston  theory  type  analysis  which  has  been  thoroughly  evaluated  : 
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in  Reference  (2).  The  question  then  arises: 


how  can  ground  vibration 


testing  on  a  "two-dimensional"  wing  lead  to  an  estimate  of  flutter  speed? 


First  place  two  concentrated  damping  forces  on  the  wing  which  produce 
the  coefficients  and  C0  given  by  piston  theory.  Then  apply  a  shaker 
force  at  the  plate  mid-chord.  What  relationship  does  that  force  have  to 
an  aerodynamic  force  Nef?  What  data  are  needed  to  compute  flutter  speed? 


* 


The  two  forces  may  be  related  by  thinking  of  the  following  experi¬ 
ment.  The  sinusoidal  ground  test  force  produces  a  sinusoidal  pitching 
angle  (Fig.  2).  It  varies  both  in  magnitude  and  phase  relative  to  the 
driving  force  as  the  frequency  is  changed.  At  some  special  frequency 
(Test  #1,  Fig.  2),  the  pitching  angle  Oc  may  be  180  degrees  out  of  phase 
with  a  downward  shaking  force.  Next,  allow  very  low  density  air  to  flow 
over  the  airfoil  at  a  chosen  supersonic  Mach  number  (Test  #2).  At  the 
instant  when  the  shaker  force  is  downward,  the  airfoil  is  pitching  nose 
down;  and  the  small  airforce  which  is  generated  is  downward,  exactly  in 
phase  with  the  shaker  force  at  the  same  point  where  the  shaker  force  is 
applied.  Under  these  circumstances,  the  sum  of  the  shaker  and  aerodynamic 
forces  increases  the  angular  oscillatior..  Reduce  the  test  shaker  force 
such  that  the  sum  of  the  shaker  and  aerodynamic  forces  and  the  resultant 
angle  of  oscillation  remain  the  same  as  before.  For  Test  #3  allow  the 
density  of  the  airflow  to  be  increased  to  a  critical  level,  where  it  is 
observed  that  the  shaker  force  is  no  longer  needed  to  maintain  the  oscilla¬ 
tion.  A  condition  of  flutter  has  been  reached. 


i  . 

9 
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The  magnitude  of  the  aerodynamic  force  is  thus  equal  to  the  shaker 
force  which  it  fully  replaces  in  causing  the  oscillation.  It  follows 
that  the  present  aerodynamic  force  per  unit  angle  of  airfoil  pitching  is 
simply  the  shaker  force  per  unit  angle.  It  is  now  seen  that  this  informa¬ 
tion  is  what  is  needed  to  make  a  flutter  speed  or  Mach  number  estimate. 

The  flutter  Mach  number  at  the  chosen  altitude  will  produce  an  aerodynamic 
force  per  unit  angle  equal  to  the  shaker  force  per  unit  angle  derived  from 
direct  ground  vibration  testing  at  zero  airspeed. 

It  is  desirable  to  underline  the  significance  in  the  experiment  of  the 
requirement  that  the  chosen  frequency  should  be  one  in  which  the  angle  is 
180  degrees  out  of  phase  with  the  test  shaker  force.  If  any  other  frequency 


Tfc*  John*  Wopfciru  Unh*f*lty 
AMM.IIO  PKYtlCf  LAM«AT«*r 
Silvor  Sp<ii^,  Maryland 


12. 


had  been  chosen,  the  phase  of  the  angle,  and  the  phase  of  the  aerodynamic 
force  due  to  the  airflow,  would  not  allow  for  the  direct  replacement  of  the 
portion  of  the  test  shaker  force  by  the  aerodynamic  force. 

The  two-dimensional  flutter  Mach  number  is,  according  to  the  fore¬ 
going  argument,  one  which  satisfies  the  following  equation. 


(15) 
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4.  Flutter  Analysis  -  Three-dimensional  System 

Both  the  theoretical  and  the  experimental  discussions  have  been  based 
on  a  two-dimensional  system.  It  is  useful  to  extend  these  concepts  to  a 
three-dimensional  system,  such  as  a  clamped  cantilever  plate. 

When  the  motion  of  the  plate  is  idealized  in  terms  of  two  degrees  of 
freedom,  primary  bending,  and  primary  torsion,  the  extension  of  the  con¬ 
cept  can  proceed  conveniently.  Lagrange's  equations  are  needed  to  define 
the  plate  motion..  In  general: 


(16) 
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where  T  is  the  kinetic  energy  of  the  system,  Vis  the  potential  energy, 
and  is  the  generalized  force  in  the  ith  degree  of  freedom. 

The  displacement  of  the  flexural  axis  of  the  plate,  namely  the  mid¬ 
chord  is  approximated  by  one  degree  of  freedom. 
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The  angular  displacement  of  the  plate  is  approximated  by  another 
degree  of  freedom. 


(18) 


a  ■  D0  q  a  »  where 


is  the  primary  torsional  mode  shape,  and  q^,  is  the  unknown  time-varying 
function  in  this  degree  of  freedom. 


The  kinetic  energy  of  the  system  is: 
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The  potential  energy  of  the  system  is: 


The  variational  work  per'  led  by  an  applied  shaker  force  P  at 
the  mid-chord  and  by  the  forward  and  aft  dampers  is: 
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When  the  shaker  force  is  absent,  and  the  surface  is  flying  at  some  super¬ 
sonic  Mach  number,  the  variational  work  performed  by  the  simplified  aero¬ 
dynamic  forces  is: 
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The  application  of  Lagrange's  Equations  yields  in  the  case  of  forced 
sinusoidal  vibration: 
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Where: 


tt  ^  is  the  uncoupled  spanwise  bending  frequency. 

»  is  the  uncoupled  spanwise  torsional  frequency. 
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In  the  case  where  the  surface  is  subjected  to  airflow  at  some 
supersonic  Mach  number,  the  application  of  Lagrange's  Equations  yields: 

.X 
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Using  the  notation: 

(25)  N  - 

it  may  be  observed  that  equations  (24)  lead  to  the  determinant  of 
Equation  (10).  It  follows  that  the  critical  value  of  N  needed  for 
instability  is  given  by  Equation  (14),  and  the  frequency  oa^,  is  given 
by  Equation  (12).  A  completely  theoretical  study  of  the  cantilever 
plate  may  thus  be  performed,  subjects '  to  the  basic  assumptions  regard¬ 
ing  the  motion  to  be  limited  to  that  described  by  two  degrees  of  freedom, 

5.  Ground  Vibration  Testing  -  Three-dimensional  Sysl.em 

The  relation  between  the  ground  test  shaking  force  and  the  distributed 
aerodynamic  forces  which  act  on  the  plate  flying  at  some  supersonic  Mach 
number,  may  be  established  in  a  manner  similar  to  that  employed  in  the 
case  of  the  two-dimensional  system.  Where  before  the  ground  test  input 
was  a  force,  the  appropriate  reference  input  now  is  that  of  generalized 
force.  In  performing  the  experiment,  a  useful  response  function  is  the 
angular  vibration  at  some  reference  spanwise  station.  Let  such  a  vibra¬ 
tion  be  found  at  a  certain  frequency  to  be  180  degrees  out  of  phase  with 
the  applied  force.  Then  allow  very  low  density  air  at  the  chosen  super¬ 
sonic  Mach  number  to  flow  over  the  surface.  At  the  instant  when  the  test 
shaker  force  is  down,  the  surface  is  pitching  nose  down,  and  the  distributed 
aerodynamic  forces  will  do  work  on  the  system  in  the  bending  degree  of  free¬ 
dom  according  to  the  expression  for  the  generalized  force. 
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In  order  for  the  level  of  the  vibration  to  remain  at  its 
initial  test  value,  the  test  input  force  may  now  be  reduced  to  a 
value  P  ,  such  that  Q-  is  the  same  as  Q.  . 
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The  value  of  q  is  recognized  by  its  effect  on  the  reference  angle 
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%  =  °ab  *  ^en  all°w  the  density  of  the  air  to  be  increased  far 

enough  so  that  the  presence  of  P  is  unnecessary  to  the  maintainance 
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The  net  generalized  force  acting  is  then  composed  of  the  sum  of 
the  shaker  and  the  aerodynamic  inputs: 


^4 
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(28)  Q.  -  -  N  q 
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This  value  of  has  remained  the  same  throughout  the  discussion  and  has 
the  value  which  was  given  it  by  the  initial  test  force.  Consequently: 


(29) 
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The  critical  value  of  N,  with  q„  deduced  from  o  is  therefore: 
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It  is  then  possible  to  define  a  critical  ratio  of  the  driving 
shaker  force  to  the  angle  of  vibration  at  the  chosen  angular  pickup 
station. 
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Using  the  classical  two-dimensional  lift-curve  slope  for  CT  ,  the 

ull 

critical  atmospheric  pressure  which  cannot  be  exceeded  for  any  super¬ 
sonic  Mach  number  is: 


(32)  p - ± - ; 

.7  CT„  M' 

La 

It  will  be  noted  that  the  grouping  of  the  terms  shows  an  aerodynamic 
factor,  a  structural  mode  shape  factor,  and  a  dynamic,  response  factor. 

The  first  has  to  be  obtained  directly  from  theory  or  from  wind-tunnel 
data,  while  the  second  and  thii*d  terms  can  be  obtained  from  theory,  from 
open  loop  ground  vibration  testing,  or  from  closed  loop  one-channel 
simulator  testing. 

6.  Flutter  Simulation 

The  open-loop  type  of  ground  vibration  test  required  that  the 
angular  response  be  monitored  by  the  test  engineer  to  find  the  frequency 
at  which  the  angle  was  180  degrees  out  of  phase  with  the  driving  force. 

To  bypass  this  particular  part  of  the  problem,  a  closed  loop  system  may 
be  set  up  which  in  effect  causes  a  shaker  force  to  be  generated  when  the 
plate  pitches,  analagous  to  the  aerodynamic  force  which  is  produced  by 
an  angular  motion  of  the  surface  in  flight.  At  some  particular  gain 
setting  of  the  closed  loop,  an  oscillation  is  initiated  at  the  appropriate 
frequency,  analogous  to  a  flutter  condition  which  would  occur  in  flight 
if  the  air  density  and  Mach  number  were  such  as  to  produce  a  certain  lift 
per  unit  angular  motion.  The  system  gain  in  the  closed  loop  case  is 
given  directly  by  the  factor  P  in  pounds  of  driving  force  per  unit 

reference  angle,  and  the  atmospheric  pressure  versus  Mach  number  which 
defines  the  flutter  boundary  would  again  be  found  from  Equation  (32). 

The  closed  loop  result  should  be  the  same  as  the  open  loop.  But  there 
are  some  differences  in  the  process  which  would  be  significant  if  more 
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than  one  shaker  were  to  be  used.  The  components  needed  for  the  closed 
loop  system  are  required  to  be  free  of  phase  lags  over  a  wide  frequency 
band.  The  transducer  gain  should  be  constant  over  a  wide  band.  Many 
components  suitable  for  getting  open-loop  data  are  inadequate  for  closed- 
loop  studies.  But  the  chief  advantages  of  the  closed-loop  simulation  would 
appear  when  more  than  one  shaker  would  be  used.  As  the  number  of  shakers 
are  increased,  the  system  converges  to  a  closer  and  closer  simulation  of 
the  distributed  airforces  which  are  produced  in  flight. 

In  this  report,  the  emphasis  on  the  one-channel  has  been  placed 
because  it  poses  severe  requirements  on  the  components,  and  serves  to 
show  their  successful  performance  under  live  conditions. 

7.  Flutter  Simulation  for  Plate  A 

A  study  has  been  performed  on  a  30”  :<  12"  x  1/8"  cantilever  plate 
described  by  the  mode  shape,  stiffness  and  inertia  data  of  Table  2.  Using 
the  data  in  the  formula  of  Equation  (8),  the  following  theoretical  result 
is  obtained  for  the  atmospheric  pressure  versus  Mach  number  above  which 
the  plate  would  flutter. 


(33) 
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An  experimental  result  was  obtained  for  a  shaker  installed  at  the 
80%  span  at  the  mid-chord.  During  self-excitation,  the  shaker  force 
was  .265  lb.,  and  the  angle  of  vibration  was  .00128  radians.  The  P  / 

fl  fl 

of  Lj7  lb/radian  was  used  in  Equation  (32), 


(34) 
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It  is  seen  that  the  one-channel  simulation  yields  a  very  unconservative 
estimate,  in  that  the  plate  is  supposed  to  be  safe  at  an  atmospheric  pres¬ 
sure  25%  greater  from  the  simulator  results  as  compared  co  a  theory  which 
incorporates  inertia  and  stiffness  data  based  upon  the  theoretical  mode 
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shapes  and  measured  frequencies.  A  number  cf  possible  explanations 
for  this  discrepancy  involve  the  questions  of: 

(1)  Using  enough  mode  shapes  in  the  theory  to  describe  the  plate 
as  it  is  driven  by  a  force  at  one  point. 

(2)  The  shaker  force  calibration. 

(3)  The  angular  pickup  calibration. 

It  is  not  believed  that  items  (2)  and  (3)  were  performed  as  accurately 
as  would  be  needed.  These  techniques  were  poorly  developed  at  the  time. 

The  shaker  which  was  used  in  the  test  was  modified  from  the  original  Goodman 
design  to  free  the  coil  from  a  20  lb/in.  spring  restraint.  Such  a  spring  was 
much  too  great  relative  to  the  spring  of  the  primary  bending  mode  at  the 

O 

shaking  point,  4. 13/ .88^  or  5.31  lb/in.  When  the  shaker  was  thus  modified, 
the  annulus  of  the  magnet  had  to  be  enlarged  to  permit  small  coil  mal¬ 
alignments.  The  force  per  amper«.  factor  then  became  a  variable  which 
required  calibration  for  each  test.  The  shaker  rod  was  preloaded  against 
the  plate  and  the  magnitude  of  the  preload  determined  by  loading  the  rod 
up  to  breaking  contact.  The  shaker  current  needed  to  cause  chatter  was 
then  measured,  to  provide  the  force  per  ampere  calibration  factor. 

The  angular  measuring  system  comprised  two  linear  differential  trans¬ 
formers  placed  at  the  leading  and  trailing  edges  of  the  plate.  A  transis¬ 
torized  phase  comparator  and  10  KC  excitation  unit  served  to  apply  the 
carrier  to  the  differential  transformer,  and  to  rectify  and  demodulate  the 
output  signal.  The  output  signal  was  fed  into  a  transistorized  power 
amplifier  which  delivered  a  current  into  the  shaker  in  phase  with  the 
difference  of  the  signals  from  the  two  pickups.  In  this  investigation 
one  of  the  problems  which  turned  up  and  which  continued  to  be  troublesome 
was  that  of  matching  the  gains  of  the  two  pickups.  It  was  soon  evident 
that  matching  the  gains  for  one  pair  of  initial  positions  of  the  pickups 
failed  for  other  pairs  of  positions,  and  closer  examination  of  the  pickup 
calibration  curves  revealed  slight  non-linearities  which  accounted  for 
the  deficiencies  in  the  system.  A  certain  amount  of  carrier  signal  could 
not  be  filtered  out  of  the  final  pickup  output,  and  it  was  evident  that 
further  amplification  of  the  output  was  impractical  for  use  in  a  higher 
gain  system.  The  phase  lag  in  the  0-200  cps  frequency  band  was  low,  but 
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any  increase  in  the  filtering  was  seen  to  increase  the  phase  lag,  and 
prevent  the  simulation  of  the  no-lag  lift  force  as  a  function  of  the 
angle  of  attack  of  the  plate. 

On  the  basis  of  the  experience  provided  in  the  flutter  simulation  of 
Plate  A  it  was  concluded  that  a  better  angular  pickup  was  needed  and  that 
it  would  be  desirable  to  use  an  unmodified  Goodman  shaker  for  those  cases 
in  which  stiff  plates  were  to  be  studied. 

8.  Flutter  Simulation  for  Plate  B 

A  number  of  plates  and  root  conditions  were  involved  in  some  further 
testing,  but  the  test  assembly  shown  in  Figure  10  was  subjected  to  the 
longest  series  of  tests  for  the  development  of  suitable  angular  pickups, 
preamplifiers  and  power  amplifiers.  These  units  were  finally  all  revised 
to  meet  the  exacting  requirements  for  the  production  of  stable  linear  out¬ 
put  force  versus  input  plate  angular  motion. 

The  detection  of  angular  displacement  provided  a  particularly  diffi¬ 
cult  problem.  It  was  natural  at  first  to  think  of  two  displacement  pickups, 
as  used  in  Plate  A.  When  difficulties  were  encountered  from  non-linearities, 
the  next  approach  was  to  devise  a  mechanical  linkage  which  vrould  subtract 
the  displacements  of  two  points  on  the  plate  spaced  a  certain  chordwise 
distance  apart,  to  produce  a  linear  displacement  proportional  to  the  mean 
angle.  The  mean  displacement  was  at  first  sensed  by  a  linear  differential 
transformer.  For  a  number  of  reasons  this  arrangement  was  not  satisfactory. 
Among  the  evils  was  the  presence  of  60  cps  noise  which  precluded  sufficient 
amplification  for  the  necessary  self-excitation  of  a  stiff  surface  with  a 
high  supersonic  flutter  speed. 

It  was  then  considered  that  the  ideal  low  noise  -  high  gain  pickup 
might  be  found  in  the  phonograph  realm,  where  such  characteristics  were 
especially  prized.  A  crystal  pickup  was  tried,  but  the  phase  characteristics 
were  poor  and  contact  of  the  needle  on  the  subtraction  linkage  was  hard  to 
maintain.  At  this  point  a  Weathers  capacitance  pickup  was  investigated. 

Its  phase  and  noise  characteristics  were  good,  but  it  was  non-linear;  so 
it  could  be  used  only  over  a  small  distance.  Precise  measures  of  the  non¬ 
linearity  were  not  obtained,  but  the  calibration  of  the  pickup  had  to  be 
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performed  frequently,  and  if  the  gap  between  the  face  of  the  pickup  and 
the  moving  surface  changed  by  .010  inches,  the  gain  changed  by  10%. 

Different  linkages  vere  employed  vi t h  this  pickup.  One  of  the  -ore 
successful  vas  a  horizontal  panel  supported  on  three  corner  feet.  Figure 
6(a).  Two  of  the  feet  ((1)  and  (2))  rested  on  two  chord vise  ptints  on 
the  specirsen  plate.  The  third  foot  (3)  rested  on  an  external  fixed  base, 
while  at  the  fourth  comer  a  Weathers  pickup  vas  placed  directly  beneath 
the  pickup.  The  pickup  vas  positioned  such  that  when  the  tvo  feet  on  the 
plate  vere  moved  equally  the  output  voltage  vas  zero.  An  output  voltage 
vould  appear,  then,  only  vhen  the  tvo  plate  feet  moved  dif ferentially. 
i.e. ,  vhen  the  plate  took  on  an  angular  displacement.  Along  vith  this  unit, 
a  voltage  preamplifier  vas  needed,  and  vhile  a  D. C.  device  vould  have  been 
desirable,  it  va3  act  feasible  vith  the  shifty  output  from  the  Weathers 
oscillator.  A  low  frequency  A. C.  preamplifier  vas  feasible  and  vas 
employed  along  vith  a  D. C.  transistorized  power  amplifier  to  produce  a 
snaker  current  in  phase  vith  the  plate  angular  displacement  over  a  vide 
frequency  band. 

Another  fairly  successful  form  of  an  angular  transducer  vas  a  vertical 
post  mounted  on  the  plate  vith  a  vertical  plane  surface.  The  Weathers 
pickup  (5) ,  Fig.  (6b)  placed  to  face  the  post  could  see  it  move  fore  and 
aft,  but  not  vertically,  and  .hus  sensed  angle.  With  the  three-footed 
panel,  questions  of  panel  modes,  faithfulness  of  feet  contact ,  friction, 
and  fixed  base  modes  all  appeared  at  one  time  or  another  and  periodically 
came  clanking  across  the  scene.  In  the  case  of  the  vertical  post-plane 
arrangement,  the  primary  source  of  difficulty  was  at  least  redjced  to  the 
modes  of  that  one  structure^vith  no  ribbing  friction  present  to  produce 
unwanted  phase  lags,  3ut  the  fundamental  non-linearity  of  the  pickup  itself 
was  always  present. 

Another  type  of  capacitance  pickup  was  considered  briefly  becaus. 
it  was  much  more  linear.  It  was  recognized,  however,  that  a  mechanical 
structure  of  some  kind  would  be  necessary,  and  would  once  nv're  provide 
modes  to  confuse  the  end  result.  For  a  long  time,  several  API  engineers, 
namely  W.  Tynan  and  R.  Hires,  had  been  suggesting  the  use  of  optical 
techniques  to  sense  the  angular  motion.  Considerai  I  r.  of  such  technique-, 
has  been  rewarding  in  closer  attainment  of  a  low  ncise^zero  phase  lag,  high 
gain,  modeless  detection  system. 
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The  system  now  comprises  as  a  light  source  a  light  beam  pointer 
of  the  type  used  in  news  meeting  presentations , in  an  arrangement  snown 
in  Figures  (6c)  and  (10).  The  beam  is  reflected  from  2  mirror  placed 
vertically  on  the  test  plate,  and  an  image  is  projected  upon  a  two¬ 
photocell  screen  some  30  inches  away.  The  cells  are  connected  in 
parallel,  with  the  positive  and  negative  terminals  at  one  end  grounded, 
and  with  the  signal  generated  at  the  junction  of  the  other  terminals* 

The  signal  is  applied  t  -.Cz  f-r-; amplifier ,  Figure  7.  The  signal  is 
zero  when  both  o'  _ne  cells  are  illuminated  equally.  The  output  of  the 
D. C.  preamplifier  is  nulled  for  the  particular  pair  of  6  volt  batteries 
used  as  a  power  source.  The  gain  control  then  does  not  change  the  output 
null.  The  signal  is  raised  to  .5  volts  to  enter  the  D. C.  power  amplifier 
Figure  8.  This  particular  circuit  is  designed  such  that  a  steady  current 
cf  .5  amperes  through  the  shaker  is  provided  by  the  input  reference  signal 
level  oi  .5  volts  from  the  preamplifier.  A  complete  D. C.  system  is  thus 
attained. 

A  complete  test  set-up  comprised  the  15.5”  x  10”  x  .25”  cantilever 
Plate  B,  clamped  to  a  heavy  aluminum  block,  and  the  one  channel  simula¬ 
tion  system.  The  simulation  system  comprised  the  electro-optical  angular 
pickup,  the  associated  preamplifier ,  the  power  amplifier,  and  a  small  3 
watt,  one  pound  force  shaker.  The  shaker  was  bolted  to  the  plate  support 
block  12"  from  the  root.  The  shaker  rod  was  then  attached  to  Plate  B  at 
the  mid-chord.  A  . 5  lb  mass  was  cemented  to  the  tip  trailing  edge  of 
the  plate,  to  provide  mass  unbalance.  The  mirror  for  the  angular  pickup 
was  attached  to  the  mid-chord  3/A"  outboard  of  the  shaker  rod  point. 

Token  damping  forces  were  supplied  by  two  1/8"  screws  located  3"  ahead 
and  3"  behind  the  mid-chord  1/2 !  in  from  the  tip.  The  ends  of  the  screws 
moved  in  cups  of  SAE  90  oil  to  produce  damping. 

The  operational  sequence  was  as  follows.  With  the  loop  gain  set 
by  the  preamp  control,  the  plate  was  tapped  lightly.  The  vibrations 
damped  out;  the  gain  was  increased,  and  the  tapping  repeated.  Eventually 
at  a  certain  gain  position,  the  plate  continued  to  oscillate.  At  larger 
gains,  the  level  of  vibration  was  such  as  l-  exceed  the  1/2  ampere  limita¬ 
tion  of  the  power  amplifier.  The  peaks  of  the  wave  were  clipped,  and  the 
system  was  in  a  non-'lir.ear  state.  At  this  point  the  oscillation  had  to 
be  stopped  so  t’.o.  the  sequence  of  events  could  be  repeated  for  the  study 
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of  the  linear  system,  (Only  cursory  examination  has  been  made  so  far 
of  the  tolerance  in  gain  position  for  which  the  oscillation  will  be 
maintained  in  the  linear  regime  of  the  feedback  loop.  The  tolerance 
is  apparently  small.) 

The  next  problem  was  to  establish  the  gain  of  the  loop  which  pro¬ 
duced  the  oscillation.  A  direct  approach  was  to  move  the  photo-cell 
pair  a  known  distance  6  and  record  the  change  in  the  voltage  E^  across 
a  one  ohm  series  resistor  in  the  voice  coil  line.  Using  the  perpen¬ 

dicular  distance  £  from  the  mirror  to  the  photocell  and  the  force  per 
P 

unit  current  the  force  per  unit  angle  was  determined  from: 


(35) 


2  X 


E 

_£ 

6 


Several  difficulties  were  evident:  first,  the  movement  of  the  cell 
had  to  be  small,  so  that  the  voltage  shift  would  not  exceed  the  ±  .5  volt 
limit  about  the  neutral  .5  volt  point.  The  movement  was  produced  by  a 
cell  support  slide  pushed  by  the  end  of  a  micrometer,  and  there  was  some 
question  about  the  fidelity  of  such  an  arrangement.  The  second  problem 
arose  from  a  shifty  voltage  output  from  the  system  in  the  high  gain  posi¬ 
tion.  To  solve  these  problems,  the  dynamic  gain  was  established  from  forced 
vibration  of  the  plate,  as  a  function  of  gain  dial  reading.  The  modified 
calibration  process  was  to  determine  the.  static  reference  low  gain.  A 
large  movement  of  the  photocell  was  necessary  to  produce  output  voltage 
shift  of  .5  volts.  Further  to  define  the  input  movement^  a  dial  gage 
reading  to  increments  of  1/10000  in.  was  placed  to  measure  S  .  The  out¬ 
put  voltage  signal  which  fluctuated  ±  .02  volts  about  the  mean  value  when 
the  gain  setting  was  high,  was  barely  perceptible  when  the  gain  was  low. 

Thus  it  was  felt  that  the  input  motion  and  the  output  signal  were  fairly 
well  defined  for  the  calibration  process. 

In  a  given  problem,  the  loop  gain  for  plate  oscillation  was  found 
by  recording  the  gain  dial  reading.  Then  the  low  gain  calibration  was 
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performed  to  find  (E/6)c>  and  the  preamp  characteristic  curve  of 
gain  ratio  versus  dial  reading  was  used  to  find  the  gain 

ratio  for  the  given  test.  The  necessary  force  per  unit  angle  was 
finally  found  from:  f- 


Equation  (32)  was  then  used  to  define  the  atmospheric  pressure  — 

Mach  number  flutter  boundary  and  the  altitudo-Mach  number  boundary 

depicted  in  Figure  9.  Calculations  for  this  purpose  are  shown  in 

4 

Table  5.  If  the  lift  curve  slope  is  given  by  — ■  —  ,  then  the 

conclusion  of  this  study  is  that  the  plate  would  flutter  at  a  Mach 
number  of  2.9  at  sea  level. 


9.  Qualifications 

At  this  point  it  would  be  advisable  to  list  the  qualifications  on 
any  result  from  the  one  channel  simulation. 

1.  The  flutter  speed  depends  upon  the  mode  shapes  and  D^.  The  shapes 
used  in  this  study  are  simply  those  for  a  uniform  cantilever.  They 
could  be  improved  by  computations  for  this  specific  mass  distribution, 
'"hey  could  also  be  measured. 

2.  The  calculation  of  the  speed  implies  that  flutter  is  taking  place  in 
a  shape  composed  primarily  of  simple  bending  and  twisting.  If  the 
second  bending  mode  played  a  significant  part,  then  the  one-channel 
representation  would  fail. 

3.  The  tip  airforces  for  this  three-dimensional  system  have  centers  of 
pressure  ahead  of  the  elastic  axis.  Again,  this  effect  is  hard  to 
incorporate  into  one  channel.  The  effect  of  forward  centers  of 
pressure  could  be  assessed  by  moving  the  shaker  further  forward. 

In  special  cases,  to  produce  flutter  speeds  on  actual  missile  wings, 
one  channel  might  be  sufficient.  But  to  produce  reasonably  good  results, 
it  is  the  writer's  opinion  that  at  least  four  channels  should  be  used  on 
surfaces  where  chordwise  bending  is  expected  to  be  unimportant.  A  crude 
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appraisal  of  the  significance  of  chordwise  bending  night  then  be  achieved 
by  the  use  of  four  more  channels* 

Ultimately  the  simulation  of  flutter  should  be  helpful  in  checking 
flutter  speeds  by  a  technique  which  is  independent  of  any  structural 
vibration  analysis  of  the  giv  n  wing.  Speed  predictions  need  not  be 
based  upon  legislated  correlations  between  theoretical  and  experimental 
mode  shapes  and  frequencies,  for  those  cases  where  frequencies  are  close 
together  and  non-linear  characteristics  abound.  The  effects  of  different 
theoretical  airforce  coefficients  can  be  appraised,  and  t^e  minimum  flutter 
speed  may  then  be  estimated. 
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TABLE  3 


Mode  Shape  Data  for  Plate  B 
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TABLE  4 


Plate  B 

Theoretical  Inertia  and  Stiffness  Data 
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These  have  been  determined  such  that  the  coupled  frequencies  agree  with  the 
measured  values  of  28  and  89.5  cps  for  the  "bending"  and  "torsion"  modes  respec¬ 
tively. 
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TABLE  5 

Flutter  Boundary  Data.  Plate  B 
Atmospheric  Pressure  Vs.  Mach  No 
Altitude  Vs,  Mach  No. 
Reference  Figure  3 


Ref.  Equation  32 
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Ref.  Equation  36 
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Planform  -  Plate  B 
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Flutter  Simulator 
Block  Diagram  -  One  Channel 
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A  -  Angular  Vibration  Transducer 
B  -  Preamplifier 
C  -  Power  Amplifier 
D  -  Shaker 
e  -  Forward  Damper 
f  -  Rear  Damper 
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FIGURE  9 

Plate  B  Flutter  Boundary 
Atmospheric  Pressure  vs.  Mach  No. 
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